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Abstract. We study the dynamics of fronts when both inertial effects and external fluctuations are taken 
into account. Stochastic fluctuations are introduced as multiplicative noise arising from a control param- 
eter of the system. Contrary to the non-inertial (overdamped) case, we find that important features of 
the system, such as the velocity selection picture, are not modified by the noise. We then compute the 
overdamped limit of the underdamped dynamics in a more careful way, finding that it does not exhibit 
any effect of noise either. Our result poses the question as to whether or not external noise sources can be 
measured in physical systems of this kind. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise and Brownian motion - 05.45.-a Nonlin- 
ear dynamics and nonlinear dynamical systems - 47.54. +r Pattern selection; pattern formation - 47.20. Ky 
Nonlinearity (including bifurcation theory) 



1 Introduction 

Front propagation is being the subject of very active re- 
search in the last few years: Indeed, the problem of the se- 
lection of the front velocity is a paradigm of the dynamical 
selection mechanisms arising in a large number of physi- 
cal, chemical and biological systems with a certain kind 
of instability (see JjJ and references therein). One of the 
important questions into which interest has been focused 
is how the deterministic front scenario is modified by the 
presence of noise. In this context, the effect of stochastic 
fluctuations on front dynamics and the modification of its 
deterministic features have been considered by several au- 
thors. A detailed summary of those results, which were 
mostly devoted to the changes of the front velocity and 
the spreading of the front position, can be found in . 

The present work addresses a related problem which, 
on the other hand, arises naturally from the above line 
of research: Is the influence of (external) noise on front 
propagation the same if inertial effects (as far as we know, 
neglected in previous work) are taken into account? As 
is well known, including inertia leads to a description in 
terms of a damped, hyperbolic partial differential equa- 
tion. It is important to note that first, a model like this 
arises when one considers a more realistic jump process for 
the individuals whose probability density is described by 
the partial differential equation; and, second, that hyper- 
bolic equations of this kind describe many actual physical 
phenomena, such as, e.g., population dynamics, nonlinear 



transmission lines, cell motion, branching random walks, 
dynamics of ferroelectric domains, and others 

The role of inertia in the scenario of (deterministic) 
frontpropagation has been studied recently JTo| (see also 
[H@j|§|j an d specially JTl| for a detailed study of the under- 
damped dynamics restricted to the linear regime). In this 
paper, it was proven that the different dynamical regimes 
of front propagation using deterministic parabolic models 
do not change, i.e., the values of the control parameter 
separating the different regimes do not depend on the in- 
ertia parameter ("mass"). However, it was also shown in 
[ fl0| that the values of the front velocities corresponding 
to every interval of the control parameter and the spatial 
shape of the propagating fronts do depend on the inertia 
parameter. 

As stated above, our explicit objective will now be the 
study of the role of multiplicative noise and its compar- 
ison to the results in |lj,[l3|], in order to understand the 
interplay of inertia and external stochastic perturbations. 
Accordingly, we undertake the study of a hyperbolic par- 
tial diferential equation with a multiplicative noise term, 
used to model front dynamics subjected to both inertia 
and external fluctuations. Opposite to the deterministic 
case, in which for very small mass or inertia a naive adia- 
batic elimination procedure (i.e., leaving the second time 
derivative term out) gives a parabolic equation which de- 
scribes very accurately the front dynamics JTo| ], we will 
prove that this is not at all the case when noise is present. 
As a matter of fact, the starting hyperbolic equation with 
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noise in the Stratonovic interpretation transforms, upon 
adiabatic elimination, into a parabolic equation with an 
extra term coming from the noise. Furthermore, the so 
obtained equation turns out to be equivalent to the usual 
parabolic equation if interpreted in the ltd sense. In ad- 
dition, we include numerical simulation results confirming 
this unexpected result. We report on these results along 
the following scheme: We begin by presenting our model 
and by briefly summarizing what is known about the noise 
influence in the overdamped case. Subsequently, we con- 
cern ourselves with the study of the externally perturbed 
(stochastic) case. We conclude by summarizing our main 
findings and discussing their implications. 

2 Model definition and notations 

Let us begin by introducing the purely deterministic prob- 
lem. Generally speaking, the situation which we are inter- 
ested in is generally modeled by the hyperbolic equation 

(j>u + OL(j) t = Dcji^ + r" 1 f^.a) (1) 

where a is the friction (dissipation), D is the diffusion 
coefficient, r is the characteristic time of the reaction term 
and a is the external control parameter of the nonlinear 
reaction term f((f>, a). The first step will be the reduction 
of the number of parameters by introducing the change of 
variables t — > rat, x — > \J tDx; our initial model, Eq. (|l|), 
reduces then to 

ecfrtt + 4>t = 4>xx + f(<P, a), (2) 

where a new parameter (the "mass"), e = (tot) , ap- 
pears. We note that the information regarding both the 
characteristic reaction time and the dissipation coefficient 
is contained now in e. With this new notation, the parabolic 
or overdamped limit is obtained by letting e->0||, which 
leads to (note that it is a singular limit) 

<th = <l>xx + f(<t>,a)- (3) 

We will refer to this procedure along the paper as naive 
adiabatic elimination. 

For the sake of definiteness, we take as a representative 
example of nonlinear reaction term 

/(^a)=0(l-0)(o + 0. (4) 

Such a term can be obtained from a local (bistable) po- 
tential, f(4>,a) — —V'(4>), with 

n« = -f^-^V+^ 4 - (5) 

It is then straightforward to show that the steady states 
are, <f>i = 0, <p2 = 1 and (f> 3 = —a. We are interested in 
those solutions which are front-like (kinks) conecting the 
(unstable if a > and metastable otherwise) state <j)\ = 
with the globally stable state 4> 2 = 1. Consequently, we 
supplement Eqs. (||) and (||) with boundary conditions 
<j>(—ao,t) = (f>2, </>(oo,t) = <j>\. 



Let us now move on to the stochastic version of the 
problem. As is well known, thermal (additive) noise is not 
expected to be relevant in actual, experimental contexts, 
as its ratio to other terms in the governing equations can 
be estimated to be 10 -9 pj. However, in addition to ther- 
mal noise we must expect [Q multiplicative noise sources 
arising from fluctuating control parameters 0]. Exam- 
ples of this case are recent experiments on the Belousov- 
Zabotinsky reaction in a light-sensitive medium 
[TtJ . The fluctuating light intensity enters as a multiplica- 
tive noise in the theoretical modelization of this chemi- 
cal reaction. According to this, noise is introduced in the 
system described so far through the parameter a, which 
fluctuates according to 

a —>■ a + £(x, t). (6) 

The noise £ is gaussian, with zero mean and correlation 
given by 

t)£(x',t')) = 2 a 2 C(x - x')6(t - t'), (7) 

with C(x — x') being the spatial correlation function, nor- 
malized by imposin g [ C( x)dx = 1. The noisy parabolic 
case considered in [|12|,[ll| corresponds to the following 
stochastic partial differential equation: 

<t>t = <f>xx + f(<f>,a)+g((l>)S(x,t), (8) 

with g{<p) — 0(1 — (f>) in case f(<fi,a) is given by Eq. (||). 
Equation (||) with noise statistical properties (0), being 
well known, will be taken as our reference scenario; its 
main features are summarized below. Nevertheless, before 
going into those results, it is important to note that, prior 
to any other consideration, it is necessary to prescribe a 
mathematical interpretation of the noise in Eq. (^) . Based 
on physical and mathematical grounds wc will follow the 
Stratonovic interpretation (see p8|,^9| for in-depth discus- 
sions of the interpretation of stochastic differential equa- 
tions). This interpretation fulfills two important proper- 
ties from a physical point of view: First, it corresponds 
to the white noise limit of a real (nonwhite) noise, and 
second, when manipulating stochastic terms, usual rules 
of calculus apply. We note that in an actual experimental 
situation the noise has necessarily a finite (non zero) char- 
acteristic time. However, this time is indeed very small 
compared with any other characteristic time of the sys- 
tem, and therefore the assumptions of white noise and 
Stratonovic interpretation seem physically well founded. 



3 Stochastic results on the overdamped 
model 

As already mentioned in the introduction, the overdamped 
limit with noise, Eq. (^), has been studied recently [p|,|T2], 
[l3| . It is worth summarizing here the main points in order 
to compare with our results. In addition, we will refer to 
the quantities introduced in these section later along the 
text. 
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In the parabolic model (g), when —1/2 < a < 1/2 
(mctastable and nonlinear regimes of the deterministic 
equation, see [20, 2l|,^2|), starting from a sufficiently lo- 
calized initial front evolves to the nonlinear solution, and 
correspondingly the selected velocity of the front is 



v n i (a) = 



2a+ 1 



V2(l-2ff2C(0)) 



(9) 



On the other hand, in the linear regime, 1/2 < a < 1, the 
velocity is given by 



t'l 



= 2y/a + a 2 C{0). 



(10) 



We note that, in the three regimes, this velocity increases 
monotonously as a function of the noise intensity; we will 
come back to this result below. 



4 Stochastic perturbation 

We now proceed to analyze the case we arc interested 
in, namely the stochastic version of the hyperbolic partial 
differential equation (En. For convenience, by introducing 
a new field variable ip, we cast it in the form 

<h = ip> 

&lh = ->l> + 4>x* + f(<t>, a) + t). (11) 

At this point, it is important to note that naive adiabatic 
elimination leads again to Eq. (|J). 

As a starting point, let us recall that in JTo| we proved 
that the velocity of the deterministic hyperbolic model 
can be obtained from the parabolic model by using the 
transformation 



V n i{e,a) 



1 



y/e + v n i (a) 



(12) 



Therefore, our first aim here is to check whether or not 
this result applies in the presence of multiplicative noise, 
i.e., whether or not we can still use the expression ( p"2"| ) 
substituting for the deterministic velocities the stochastic 
ones found in |l2|,^3). From those papers, we know that 
the important noise effects come from the fact that the 
multiplicative noise term has a non zero mean value, pro- 
vided the Stratonovic interpretation is used. Within that 
interpretation, it can be found by using only usual stochas- 
tic calculus and a little algebra that the mean value of the 
noisy term in Eq. (O) is M 



dg{4>) 6<f> , 
dcj> 5£(x',t')' 



-X, t'—t 



- 0. 
(13) 

This result allows us to conjecture that in the hyperbolic 
case the noise is not as relevant as in the parabolic case. 
In order to verify this conjecture, we carried out numer- 
ical simulations of Eq. (|ll]) (with noise white in space, 
implying C(0) = 1 /Ax with Ax the spatial discretization 
step) for different choices of the parameters, finding that 



the velocity is not affected by the perturbation even for 
very small values of e. In Fig. [l] we present some of these 
numerical results. It can be clearly seen the dramatic dif- 
ference between the effects of noise in the two models. It 
is most important to stress that the results in Fig. [l] corre- 
spond to individual realizations of the noise; of course, we 
have verified that this is the typical behavior by repeat- 
ing the simulations very many times. This means that our 
conjecture that the noise effects and in particular those 
of the front velocity are not important, is in fact true be- 
yond mean values, i.e., individual fronts do not change 
their deterministic velocity in the presence of noise. For 
the parabolic case we see that the velocity does depend on 
the intensity of the noise, being independent in the hyper- 
bolic case. This is an unexpected result but, as we have 
just seen, it is in complete agreement with our theoretical 
analysis, that leads to Eq. ([TJ]). 

Let us now discuss the reasons for this result. It is 
clear from Eq. (|l3| ) that such a null result comes from the 
fact that the response of the field to the noise is zero at 
t' = t. This is not the case in the parabolic case (|J), for 
which one finds u2, 0| 



(g(MM) = a 2 c(o)( 



dej) 5£(x>,t') 



x' =x ,t'=t 



^(0)<^#W)> 



(14) 



This non zero mean average is the responsable of the 
strong effects of noise in the parabolic model (^) which 
were studied and quantified in p^,p^[. As it turns out, 
this is not the case in the hyperbolic model (pd|); there- 
fore, it seems likely that the naive adiabatic elimination 
procedure, which led to Eq. ( g|), is not the proper way to 
take the limit e — * in Eq. (111), in view of the different 
behavior of the two models. 

In order to check this idea, we have followed an alter- 
native, non naive adiabatic elimination procedure |po| to 
see whether an equation different from (|§|) arises for the 
overdamped limit. We will outline here the main steps of 
the calculation following [^3| (see |l^] for an alternative 
presentation). For the sake of simplicity, we rewrite Eq. 
(|ll|) in a more compact form: 



1 



<Pt=i>] & = -(F(ct>)+g(<P)Z(x,t)) 



(15) 



where F((f>) — /(</>, a) + 4> X x- Formally integrating the sec- 
ond expression in ( |l5| ) and using the first one, we find an 
integro-differential equation for the variable 4>: 



e -(t-*')A 

1> = <k= I dt (F(</>)- 

h e 



g m(x,t)), (i6) 



with the initial condiction "0(0) = 0. Subsequently, formal 
integrations by parts are performed in order to obtain a 
series of terms in powers of e, which is the situation we 
are interested in (e small). From this formal expansion a 
Fokker-Planck equation, whose first order term does not 
depend on e, is obtained. The calculation is involved, but 
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Fig. 1. Front velocity versus the effective noise intensity for 
e = [circles, parabolic case (g)], e = 0.1 (squares) and e = 1 
(triangles). Filled circles correspond to the dynamics given by 
Eq. (TL7h. In all cases, a = —0.1; similar results are obtained 
for other values of a. Lines correspond to Eqs. ([]), ([jij, and 
(p^), the three straight lower lines with a = 0; among these, 
the dashed line is the theoretical prediction for Eq. 



it does not require any further physical assumptions, it 
is only (lengthy) algebra. We then skipthe details and 
refer the reader interested in them to ]23|. The final re- 
sult is that the corresponding Langevin equation in the 
Stratonovic interpretation is 



-/(<M-a 2 C(0)^.gW>)- 



■g{<f>)^x,t) (17) 



We have thus obtained a different overdamped limit, for 
which one can check that, according to Eq. (El), the mean 
value of the noisy term compensates the new term in Eq. 
(O), hence rendering the noise contribution null, as in the 
e 7? case. Therefore, Eq. (fl7| ) is the physically consistent 
overdamped limit of Eq. (|1 1|) , insofar it exhibits the same 
behavior as this last one does for any value of the "mass" 
e. 



5 Discussion and conclusions 

In this work, we have shown analytically and numerically 
that inertial effects of any magnitude suppress the external 
noise influence on the velocity of fronts. Whereas the theo- 
retical result has been conjectured by taking averages, our 
numerical simulations show that the velocity of individual 
fronts is unchanged by noise. This means that the over- 
damped (parabolic) equation usually employed to describe 
front propagation in reaction diffusion model systems is 
not simply the limit of an underdamped (hyperbolic) ver- 
sion, as in that case it is known |l2],[F| that the velocity 
of the front does depend on the noise strength. In other 
words, the naive prediction based on deterministic results, 
Eq. (|l2|), is not correct in the presence of multiplicative 



noise. We have also shown that, by means of a more in- 
volved adiabatic elimination procedure, it is possible to 
obtain an equation for the overdamped limit which does 
not show noise effects. However, this equation differs from 
the usual one by an extra term, arising from elimination, 
which exactly cancels the noise contribution to the mean 
velocity of the fronts. 

From the physical viewpoint, these results are very rel- 
evant. Indeed, we have seen that any amount of inertia 
present in the system will lead to front propagation at 
the deterministic velocity even in the presence of exter- 
nal noise. Although the calculation has been done for a 
specific model, the reasons for the vanishing of the noise 
contribution are generic and do not depend on the specific 
choice of the reaction term. In this context, it is then clear 
that, even if a system is considered overdamped, generally 
speaking there will be some degree of inertia in its dynam- 
ics. In that case, the predicted changes in the velocity [ fl2[ 
|l3f will not actually occur. Therefore, our results arc in 
fact a criterion to establish whether a system (where front 
propagation arises) is truly overdamped and correspond- 
ingly described by a parabolic equation, or in turn, it is 
an inertial system with very large damping: As we have 
seen, the response of the system to external noise would 
be fundamentally different in both cases. This result is of 
particular importance in the study of excitable media in 
noisy environments JT5|,p^7|. An additional implication of 
our findings is that, if one is interested in identifying or 
measuring possible noise effects in hyperbolic problems in 
the class of Eq. ([ll]), the front velocity is not a good ob- 
servable and it is necessary to resort to other indicators, 
such as, e.g., the fluctuations of the velocity. 

Finally, another remarkable fact is that the naive para- 
bolic limit (||) if interpreted in the Ito sense, it is stochas- 
tically equivalent to the correct limit ([Tt]) valid for the 
Stratonovic interpretation. This is but a further indica- 
tion that the overdamped limit of Eq. ( |TT| ) is problematic 
and has to be carefully performed, as in any other instance 
where, for physical reasons, multiplicative noise has to be 
considered. On the other hand, the fact that Eq. (j|) is the 
consistent overdamped limit in the Ito interpretation has 
to do as well with the known result that, in general, per- 
turbative expansions such as the one needed in the adi- 
abatic elimination procedure have many identically zero 
terms when carried out in Ito sense [fill . This result poses 
questions of a more mathematical character that would be 
interesting to address in a general framework for stochas- 
tic partial differential equations. 
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